Bose condensates at high angular momenta 
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We exploit the analogy with the Quantum Hall (QH) system to study weakly interacting bosons 
in a harmonic trap. For a ^-function interaction potential the "yrast" states with L > N(N — 1) 
are degenerate, and we show how this can be understood in terms of Haldane exclusion statistics. 
We present spectra for 4 and 8 particles obtained by numerical and algebraic methods, and demon- 
strate how a more general hard-core potential lifts the degeneracies on the yrast line. The exact 
wavefunctions for N = 4 are compared with trial states constructed from composite fermions (CF), 
and the possibility of using CF-states to study the low L region at high TV is discussed. 
PACS numbers: 03.75.Fi, 05.30. Jp, 73.40.Hm 



The close relation between high angular momentum 
states of a condensate of weakly interacting hard core 
bosons and the Quantum Hall (QH) effect was 

recently pointed out The essential observation 

is that the weak interaction limit allows for a two- 
dimensional description of the boson system in terms of 
lowest Landau level (LLL) wave functions [Q - just as for 
a QH system, and they are both described by wave func- 
tions containing powers of the Laughlin-Jastrow factor 
IXkjO^ — z i)i where Zi is the complex coordinate of the 
i th particle. As pointed out by Cooper and Wilkin ||, the 
two systems can in fact be mapped onto each other by a 
standard Leinaas-Myrheim transformation |Q , attaching 
an odd number of units of statistical flux to each particle. 
This enables us to use both intuition and techniques from 
the QH system to study rotating Bose condensates. In 
this paper we shall mainly discuss the angular momen- 
tum region N(N-l) < L < 2N(N-l), where the ground 
state corresponding to a 5-function two-body interaction 
is degenerate. We explicitly show how these degenera- 
cies can be understood via a mapping to a system of free 
anyons in the LLL, and then show how the degeneracy 
is broken by a more general short range potential con- 
taining derivatives of delta functions. We also make a 
detailed comparison between algebraically calculated ex- 
act wave functions and trial wave functions formed from 
so called compact states of composite fermions, a con- 
struction orginally due to Jain and Kawamura M. In 
particular, we will emphasize the importance of a cer- 
tain class of wavefunctions where the polynomial part is 
translationally invariant. 

Although for computational reasons we have results only 
for few particles, N — 4,6 and 8, it is clear that some 
of our results, like the degeneracy structure of the yrast 
line, hold for any N. We also believe that many features 
of the results for the CF wave functions will generalize 
to higher N . 

Since the flux attachment changes the angular momen- 
tum by mN(N — l)/2, where m is the number of fluxes 
attached, the boson - fermion mapping would apparently 
only be useful for studying angular momenta that are 
out of reach of present experiments (which are lim- 
ited to the strong interaction regime and L ~ N). How- 



ever, there are some indications that fcrmionic techniqes 
could be useful for L as low as N, i.e. for the so called 
single vortex state. Although we will return to this point 
at the end of the paper, we shall for now, without any 
further apologies, consider the theoretical problem of un- 
derstanding the region N(N - 1) < L < 2N(N - 1) of 
the yrast line. 

The simplest model for a hard-core interaction is a 
delta function potential. We thus consider a model of N 
interacting spinless bosons in a harmonic trap of strength 
lu. In the limit of weak interaction, this may be rewritten 
H as a two-dimensional lowest Landau level (LLL) prob- 
lem in the effective "magnetic" field B e ff = 2mu) with 
the Hamiltonian taking the form 



H = loL 



i<j 



(1) 



(h = 1), where L = Y]. lj — L z is the total angu- 
lar momentum. The single-particle states spanning our 
Hilbert space are rjo, i = (2 1+1 ttU)~ 1 ^ 2 z l exp(— zz/4) with 
z = \/2rnuj(x + iy). 

In Figs, g and | we show the interaction energy, in 
units of g/47r as a function of the total angular momen- 
tum L for N = 4, L < 20 and N = 8, 30 < L < 56, 
respectively. The many-body states are obtained from 
Lanczos diagonalization suitable for large and relatively 
sparse matrices [|o| . The Fock space is spanned by single- 
particle states that are characterized only by the positive 
quantum numbers I, where < I < L. (Similar exact di- 
agonalization studies have recently been reported in |J 
and pT|.) We note the following properties in the many- 
body spectra: Since increasing the angular momentum 
spreads out the particles in space, the yrast energy, i.e. 
the lowest possible interaction energy for given L, de- 
creases with increasing L. For each state in the spec- 
trum, there exists a set of "daughter states" with higher 
values of L, having exactly the same (interaction) energy 
as the original state. These daughters are simply center- 
of-mass excitations of the original states, thus having the 
same many-body correlations p2[. 
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At L = N(N — 1) there is a unique state with zero 
interaction energy corresponding to S(z\, z 2 , z^) = 1, 
while the states at higher L typically are degenerate. The 
systematics of these degeneracies can be understood by 
a mapping to anyons in the lowest Landau level. The 
essential observation is that the wave functions (||) de- 
scribe anyons in the LLL with statistics parameter a = 2 
|l3| (in general, the statistics parameter is given by the 
exponent of the Jastrow factor). It is known ]l3| , p^ , p^[ 
that anyons in the LLL obey Haldane's fractional exclu- 
sion statistics (FES) @, and following Ref. @, one can 
use this knowledge to construct the allowed many-body 
states for given N and L as angular momentum excita- 
tions of the Laughlin-like state at L = N(N—1). Accord- 
ing to the definition of FES, each particle in the system 
blocks a — 2 single-particle states, and many-particle 
states with total angular momentum L are constructed 
by occupying single-particle states, with a minimum dis- 
tance of a = 2 between each pair of occupied levels (for 
an example, see Fig. [}]). The number of allowed configu- 
rations then gives the degeneracy of the state for a given 
L. 



FIG. 1. Many-body spectra of N = 4 weakly interact- 
ing bosons in a harmonic trap for L < 20. The dashed 
line connects the yrast states 



5 




- 



30 35 40 45 50 55 60 
Angular Momentum L 

FIG. 2. As Fig. 0, but for N = 8 and 30 < L < 56, 
showing only the 100 lowest eigenvalues 




FIG. 3. FES construction of yrast states for N = 4, 
L = 12(= L ), 13 and 14. The L = 12 and L = 13 states 
are non-degenerate, whereas L = 14 has degeneracy 2. 

Table 1 shows the degeneracies of some of the states on 
the N = 4 yrast line, as obtained from the anyon map- 
ping, and they are in exact agreement with our numerical 
results. 
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For L > N(N — 1) there are zero energy states, which 
can be understood by noting that any wavefunction of 
the form 

ip(zi, z 2 , ...,z N ) = - Zj) 2 S(zi,z 2 , ...,z N ) , (2) 

where S(z±, Z2, Ziv) is a symmetric polynomial in 
the Zi'.s, has zero interaction energy. Since the fac- 
tor Ili<j( z i — z j) 2 contributes an angular momentum 
Lo = N(N — 1), states of the type (||) exist for L > L$. 



TABLE I. Degeneracy d of the lowest //-excitations above 
the Laughlin state for iV = 4. 
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This construction implicitly uses that all the eigen- 
states on the yrast line for L > N(N — 1) contain the 
Jastrow factor rii<j( z i — z j) 2 ' so tne degeneracies can 
also be found as the number of ways one can distribute 
M = L — Lq units of angular momentum among N par- 
ticles. 

For a more general hard core potential, the E = 
states on the yrast line above L = N(N — 1) are no 
longer degenerate. To demonstrate this point and study 
how the degeneracy is broken, we add a potential of the 
form 



V = Cl V 2 <5 2 (z - z') + c 2 V A 5 2 (z - z'), 



(3) 



that was originally used by Trugman and Kivelson |T^ ] 
in the context of the fractional QH effect. The term 
~ V 2 <5 2 (z — z') does not contribute to the interaction en- 
ergy for fully symmetric states, whereas the V 4 5 2 (z — z') 
term gives small corrections to the spectra in Figs. [j] 
and H (at the percent level for the parameters used in 
the inset of Fig. ^, where we show regularized forms of 
the potentials (|l|) and (||)). 

We have examined how the potential (||) splits up the 
degeneracies of the zero interaction energy yrast states, 
by exact algebraic diagonalization, using computer al- 
gebra. Here we have directly used the form (|2|), and 
systematically exploited that for a given L, all states 
corresponding to center-of-mass excitations of lower L- 
eigenstates are orthogonal to the subspace of "new" 
states. The latter subspace consists of translation in- 
variant (TI) polynomials, i.e. functions invariant under 
a simultaneous, constant shift Zi — > zi + a of all the co- 
ordinates |Q . Following Trugman and Kivelson [Q , we 
have used a basis constructed from elementary symmet- 
ric functions s n . For given N and L, the basis consists 
of all possible combinations 



\hk 2 ...k n ) = s£(z i )8%(z i )...8 k 1 ?(z N ) 



(4) 



such that 53n=i = L — Lq. Note that, for n > 2, 
we have introduced the new variables Zi = Zi — z c , with 
z c the center-of-mass coordinate z c = {^2zi)/N. The 
basis states spanning the TI subspace are identified as 
those with k\ = 0. The diagonalization is thus per- 
formed within this subspace only, which reduces the ma- 
trix dimension substantially. The resulting spectrum for 
N = A, with the coefficient c 2 in (|^) set equal to 1, is 
shown in Fig. ||. We notice the close similarity between 
Fig. [I] and f§ In both cases, the yrast line passes through 
the same number of steps, with the same step lengths, 
as L is increased by N(N — 1), to the point where the 
yrast energy becomes zero. At the point L = 2N(N — 1), 
the zero-energy yrast state is again non-degenerate and 
of the Laughlin type, i.e. Yii<j( z i ~ z j) 4 ' whereas the 
subsequent yrast states have degeneracies corresponding 
to FES with statistics parameter a — 4. Including even 
higher derivative terms in the repulsive potential would 
subsequently split up the higher regions of the yrast line. 
Finally, note that the yrast states corresponding to cusps, 



i.e. states that are followed by a "plateau" in the yrast 
line, are always in the TI subspace. 
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FIG. 4. Yrast spectrum for TV = 4 and 12 < L < 26, 
for a repulsive potential V 4 5 2 (z — z'). Translation invari- 
ant eigenstates are denoted by diamonds, whereas the 
crosses denote center-of-mass excitations. Inset: Delta 
function potential (regularized as a gaussian) and a po- 
tential including the terms (|J), with c\ — 0.023 and 
C2 = 0.00008. We see that the latter more closely re- 
sembles a hard sphere potential. 



The algebraic diagonalization used here is limited in 
practice to smaller particle numbers and angular mo- 
menta than the numerical scheme used in Figs. ^ and |^. 
However, the present approach has the advantage that 
it provides explicit, analytic expressions for the eigen- 
f unctions, in terms of symmetric polynomials. This gives 
some additional insight into the structure of the yrast 
states, and in particular the region below the single 
vortex in the case of a pure delta function interaction. 
Bertsch and Papenbrock |ll| recently proposed and nu- 
merically tested the following form for the yrast states at 
2 < L < N, 



4>L{zi) 



1 <p 2 <.. 



(Zp ± z c )[zp 2 z c )...(z. 



'Pi 



=) • (5) 



This is just the symmetric polynomial Si(zi), i.e. the 
state |0...1..0) (with the 1 in the Lth place), in the no- 
tation of (H) (without the Jastrow factor in the present 
case of a pure delta function interaction). This state is 
a basis state in the TI subspace for all 2 < L < N. In 
the cases where this is the only basis state (L = 2, 3 for 
N > 3), it is thus obvious that (||) is exact. Furthermore, 
performing the algebraic diagonalization up to L = N for 
N = 4 and N = 6, we have confirmed that even when the 
translation invariant subspace is spanned by more than 
one basis vector, the basis state (||) is always an exact 
eigenstate. 
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We now turn to a study of a class of wave functions 
that can be constructed in analogy with the so-called Jain 
states for the fractional QHE ]l7|] . The main idea is to 
map the strongly interacting LLL bosons to weakly in- 
teracting composite fcrmions by attaching an odd num- 
ber m of flux quanta to each particle. Trial wavefunc- 
tions with angular momentum L are thus constructed 
as non-interacting fermionic wavefunctions with angular 
momentum L — m/2 N(N — 1), multiplied by m Jastrow 
factors and projected onto the LLL, 

4>l=V [h{z h Zi) n^i-^rj • (6) 

Here, fp is a Slater determinant consisting of single- 
particle wave functions rj ni i(z,z) oc z L n {zz/2) where n 
is the Landau level index ( I > —n), and L n a generalized 
Laguerre polynomial. 

Originally used for the homogeneous states relevant for 
the fractional QHE, wavefunctions of the type (||) were 
later employed to describe inhomogeneous systems such 
as quantum dots [p|, ^8Ul9| ] and recently by Cooper and 
Wilkin [0] to study the bosonic yrast lines for up to 10 
particles in the case of a pure delta function interaction. 

In short, the LLL projection in (|^) amounts to the re- 
placement Zi — > 2di in the polynomial part of the wave 
function. However, there are several ways of doing this 
in practice, and we shall compare the different projec- 
tion methods when constructing trial wave functions for 
the yrast states in Figs, [j] and [|. The most straightfor- 
ward way is to replace the z:s with derivatives in the final 
polynomial, obtained after multiplying out the Slater de- 
terminant and the Jastrow factors and moving all z:s to 
the left. In practice, this method is applicable only for 
small particle numbers, when the number of derivatives 
involved is not too large. We shall refer to this method 
as "Method I" and apply it in a few special cases for 
comparison. 

Alternatively, noting that jl7) 
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where rjij = Vi( z j, z j) and Vij = Vi{zj,Zj) U. k ^ji z j~ z k) p > 
one can first absorb 2p Jastrow factors in the Slater de- 
terminant, and then project entry by entry. Since the 
wave function @ contains an odd number m of Jastrow 
factors, one finally has to compensate by multiplying the 
resulting polynomial by nj<j( z i — Zj) m ~ 2p . We shall re- 
fer to this as Method II and use it as follows; for a delta 
function interaction, where the wave function contains 
m = 1 Jastrow factor, it is appropriate to use p = 1 
(Method Ila). In the case of the V 4 (S-potential, the Jain 
construction will involve absorbing m = 3 flux quanta in 
the wave function, and we shall compare projection with 
p = 1 (Method Ila) and p = 2 (Method lib). (Note that 



Method Ila is in a sense trivial, since it relates the wave- 
functions at L and L + N(N — 1) by multiplication with 
two Jastrow factors.) 

We have already stressed the significance of the TI sub- 
space - these are the states that determine the shape of 
the yrast line. It is very appealing that there is a special 
set of the states (g) that are in this subspace, namely the 
so called compact states [Q which are characterized by 
having the n:th Landau level occupied from l n = —n to 
l n = l™ ax without any "holes" . In the context of the QH 
effect, the important property of the compact states is 
that they are homogeneous. When describing quantum 
dots using the non-interacting composite fermion model 
(NICFM), one can show that CF candidates for the cusp 
states must be compact ||, and the same line of argu- 
ments was later used for bosons From the point 
of view of the wavefunctions, the importance of the com- 
pact states is that they are in the TI subspace, and it is 
rather remarkable that for all L where a compact state 
can be constructed, the one with the lowest CF energy 
has a large overlap with the lowest exact state in the TI 
subspace. This is true, independent of whether or not the 
state is a ground state, i.e. corresponds to a cusp in the 
yrast line. In the following discussion, we shall only be 
concerned with the wavefunctions, and will not discuss 
whether or not the NICFM can give a good description 
of the energy spectrum, a que stion that was already dis- 
cussed in some detail ||) [jj| , [19[ . 

(5-function potential: In this case we have calculated 
the Jain wavefunctions ([|) corresponding to all compact 
states for N = 4, < L < 12, taking m = 1 and using 
projection method Ila (p = 1). If there are two compact 
states with the same L, we use the one with the lowest 
CF effective energy. In table || we show the overlap with 
the exact algebraic wavefunctions. A "1" indicates that 
the wavefunctions are identical. For comparison we have 
also included the overlaps with the wavefunctions cor- 
responding to projection method I, as given by Cooper 
and Wilkin || . We note that all the compact states have 
very large overlap with the exact eigenstates. That the 
CF wavef unction is exact for L — 2, 3 is a simple conse- 
quence of the TI property and that there is only one state 
in the TI subspace at these L values. That the L = 7 
state comes out exact is more surprising, and we have 
no good explanation for this. It would be interesting to 
pursue the exact diagonalization to higher N in order to 
see if there are other non-trivial CF states that are exact. 
We also note that the direct projection (method I) does 
slightly better than method II. 
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TABLE II. Overlaps between trial and exact wave yrast 
functions for iV = 4, < L < 12 and a pure delta function 
interaction, using projection methods Ila and I (see text). 
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V 4 <5-function potential: Here we calculated the Jain 
wavefunctions (|6|) corresponding to all, lowest CF-energy, 
compact states for N = 4, 12 < L < 24 taking m = 3 
and using projection methods Ila (p — 1) and lib (p = 2). 
For comparison, we also used method I for three L val- 
ues. In table III we show the overlap with the exact alge- 
braic wavefunctions. We again note that all the compact 
states (except L = 17 with projection method I) have 
very large overlap with the exact eigenstates, and that 
for L = 17, 19, the CF state is not the ground state. 
Comparing the different projection methods, we see that 
in many, but not all cases they give identical results. In 
particular, one should note that method lib does not re- 
produce the exact wavefunction for L = 14, 15 where the 
TI subspace is non-degenerate. This means that that 
projection gives a wave function that is not in the space 
of states given by (||). The same effect is even more strik- 
ing for L = 17, where method lib gives the exact wave 
function, whereas direct projection (method I) gives a 
rather poor overlap, indicating that a large component is 
not in the subspace (0). 
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TABLE III. Overlaps between trial and exact yrast wave 
functions for N — 4, 12 < L < 24, using two different pro- 
jection methods (see text). The asterisk indicates that the 
L — 17 and L = 19 Jain states (which correspond to the low- 
est TI states) are not yrast states, but lie higher in energy 
than the CM excitation of the previous L state. Note that 
projection method Ila gives zero for L — 17. 



Finally, we comment on the possibility to use CF wave- 
functions in the experimentally more relevant case of low 
L, and in particular for the single vortex state at L = N. 
Rather surprisingly, the overlaps between the CF and the 
exact wavefunction tend to get larger with increasing 
particle number, at least up to N — 10. [|| It thus seems 
worthwhile to use the CF approach to construct trial 
wavefunctions for the single vortex for general N. The 
single vortex CF state is in fact unique, if one demands 
that in addition to being compact, it should also have 
minimal CF cyclotron energy. The relevant Slater de- 
terminant is formed from the single particle states ffo,— n 
for n = N — 2, N — 1, ...0 and 7701 , and, using projection 
method I, the resulting wavefunction takes the following 
rather compact form, 



N 



L=N 



Et- 1 )" 2 " II (dk-dOllizi-Zj). (8) 

n=l k<l;k,l^n 



Although the ~ N 2 derivatives make it diffi- 
cult to evaluate this function for large N, it 
can easily be handled in integrals of the form 
J I\ t d 2 Zi exp(-l/2 J2iZiZi)f(zi)ipL=N by partial inte- 
grations. Alternatively one can use projection Method 
II, where the wavefunction becomes a determinant of lin- 
ear combinations of elementary symmetric polynomials. 
In both cases, it should be possible to compare with the 
exact wavefunction (|^) using Monte Carlo methods. 
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